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We investigate the characteristic effects of nuclear motion on attosecond transient absorption 
spectra in molecules by calculating the spectrum for different model systems. Two models of the 
hydrogen molecular ion are considered: one where the internuclear separation is fixed, and one 
where the nuclei are free to vibrate. The spectra for the fixed nuclei model are similar to atomic 
spectra reported elsewhere, while the spectra obtained in the model including nuclear motion are 
very different and dominated by extremely broad absorption features. These broad absorption 
features are analyzed and their relation to molecular dissociation investigated. The study of the 
hydrogen molecular ion validates an approach based on the Born-Oppenheimer approximation and 
a finite electronic basis. This latter approach is then used to study the three-dimensional hydrogen 
molecule including nuclear vibration. The spectrum obtained from H 2 is compared to the result of 
a fixed-nuclei calculation. In the attosecond transient absorption spectra of H 2 including nuclear 
motion we find a rich absorption structure corresponding to population of different vibrational states 
in the molecule, while the fixed-nuclei spectra again are very similar to atomic spectra. We find that 
light-induced structures at well-defined energies reported in atomic systems are also present in our 
fixed nuclei molecular spectra, but suppressed in the H 2 ’^ and H 2 spectra with moving nuclei. We 
show that the signatures of light-induced structures are closely related to the nuclear dynamics of 
the system through the shapes and relative arrangement of the Born-Oppenheimer potential energy 
curves. 

PACS numbers: 33.80.Wz, 33.20.Tp 


I. INTRODUCTION 


Recent progress in laser technology has made it possi¬ 
ble to produce isolated attosecond pulses [T], and with 
these electron dynamics can be followed on its natu¬ 
ral time scale [2]. One of the promising methods to 
follow sub-femtosecond dynamics in electronic systems 
is attosecond transient absorption spectroscopy (ATAS), 
which offers an all-optical approach to light-matter inter¬ 
actions, with the high time resolution of the attosecond 
pulse and the high energy resolution characteristic of ab¬ 
sorption spectroscopy. Since 2010, ATAS has proven it¬ 
self as a very efficient tool in studying electron dynamics 
in atoms. By following the movement of electron dis¬ 
tributions in a time integrated manner, ATAS experi¬ 
ments have been used to study a number of dynamical 
processes. Examples hereof are the observation of va¬ 
lence electron motion in krypton ions [3], autoionization 
in argon [3], the observation of AC Stark shifts of ex¬ 
cited states in helium and krypton [5], and very re¬ 
cently observation and control of the dynamics in a two- 
electron wave packet of helium [7]. Common for these ex¬ 
periments is that a near-infrared (NIR) or infrared (IR) 
pulse dresses the system, while the modification of the 
spectrum of a much shorter extreme ultraviolet (XUV) 
attosecond pulse is measured. There is some ambigu¬ 
ity in the literature about which pulse should be refer¬ 
eed to as the pump (probe) pulse in the pump-probe 
setup of ATAS. We circumvent this problem by refer¬ 
ring to the two pulses simply as the ’NIR pulse’ and the 
’XUV pulse’. A variable delay between the XUV and 


NIR pulses makes it possible to introduce an indirect 
time-resolution in time-integrated ATAS spectra, and in 
this way follow the evolution of the NIR field-dressed 
system. Theoretical investigations of ATAS have focused 
on explaining some of the typical features in the spec¬ 
tra such as Autler-Townes splittings of absorption lines 
[SI m , light-induced structures (LIS) [10] , and different 
interference phenomena [11] [12]. Often simplified few- 
level models can capture the essential physics of a pro¬ 
cess and one can benefit from the fact that the moderate 
intensity of the XUV pulse places its interaction with 
the quantum system in the perturbative regime [5]. Fi¬ 
nally an increasing number of direct calculations of the 
delay-dependent ATAS intensity spectrum has been car¬ 
ried out by solving the time-dependent Schrodinger equa¬ 
tion (TDSE) for a variety of systems. Examples hereof 
are given in Refs. Bun] where the ATAS intensity spec¬ 
tra were found by solving the TDSE for helium. In the 
latter case the single-active-electron approximation was 
used. 

Very recently ATAS experiments have been carried out 
on nitrogen m and hydrogen [l4| molecules. Neglect¬ 
ing the slow molecular rotations, the main difference be¬ 
tween atoms and di-atomic molecules is that molecular 
systems vibrate in the internuclear separation. In view of 
the emerging experimental interest in ATAS of molecules 
it is important to characterize the main effects associ¬ 
ated with the nuclear vibration. Such a characteriza¬ 
tion is the aim of the present work. First we investigate 
one-dimensional (ID) H 2 ^ by solving the TDSE exactly 
within the limitations of our two models: one where the 
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nuclear distance is fixed and one where the nuclear de¬ 
gree of freedom is included. We note in passing that 
there is a long history in strong-held physics for using 
co-linear H 2 ’** as a generic system to highlight molecular 
effects [HHni- The TDSE calculations in H 2 ^ show that 
the dissociative nature of the excited electronic states in 
the moving-nuclei model leads to major differences in the 
ATAS spectrum when compared to the hxed-nuclei case. 
The hxed-nuclei spectra are very similar to atomic spec¬ 
tra. In both cases the spectra are dominated by narrow, 
well-separated, absorption features centered around en¬ 
ergies corresponding to bound states of the system. In 
moving-nuclei spectra, on the other hand, the absorp¬ 
tion features are much broader and not necessarily well- 
separated. In the analysis of these broad moving-nuclei 
induced absorption lines we establish a method to de¬ 
termine the dissociation time of nuclear wave packets 
evolving on certain excited BO surfaces. This method 
does not depend on a pump-probe setup, and informa¬ 
tion about the dissociation time can be obtained from 
single-pulse experiments. We also hnd that the ATAS 
spectra obtained from TDSE calculations can be repro¬ 
duced accurately using an approach based on the Born- 
Oppenheimer (BO) approximation and an expansion in a 
hnite electronic basis. This approach is then used to cal¬ 
culate the ATAS spectrum of three-dimensional (3D) H 2 
including nuclear motion. Compared to H 2 ^, H 2 repre¬ 
sent a case where the lowest excited electronic states are 
binding. In the ATAS spectrum of H 2 we see a rich struc¬ 
ture associated with population transfer to vibrational 
states of excited BO surfaces. In the spectrum of H 2 we 
recognize some of the characteristics found in fixed nuclei 
H 2 and atomic spectra, however, other effects such as 
Autler-Townes splittings and LIS at well-defined energies 
are absent in the moving nuclei spectrum. Our analysis 
show that the suppression of such LIS in moving nuclei 
systems are directly related to the molecular dynamics 
of the system, and that LIS in molecules in general are 
extremely dependent on the population of the involved 
vibrational states and thereby on the BO curves of the 
involved electronic states. 

The paper is organized as follows. In Sec. |TT] we 
describe a theory for calculating the ATAS intensity 
spectrum. We use the single-system response model of 
Ref. [IH] (see also Ref. [H]). Within the approximations 
given in Sec. the ATAS intensity spectrum is readily 
determined if the dipole moment of the system is known 
for all times. In Sec. El we also establish the relation¬ 
ship between our method and other methods often used 
to calculate the ATAS spectrum [71 UHl HI]- In Sec. |Tn] 
we present calculated ATAS spectra for our two models 
of H 2 ^, and discuss the atomic-like characteristics of the 
model with fixed nuclei. A three-level model description 
capturing some of the interference features is discussed 
in Appendix]^ In the model of H 2 ^ where the nuclei are 
allowed to move, we find that the absorption lines in the 
ATAS spectrum become very broad. In Sec. |IV| we set 
up a model to explain these broad absorption lines, and 


use it to closer analyze the dissociation process of excited 
electronic states. In Sec. we present the ATAS spec¬ 
trum for 3D H 2 for fixed and moving nuclei. The fixed 
nuclei spectrum is very similar to the atomic-like, fixed 
nuclei spectrum of H 2 ^. However, as a result of pop¬ 
ulation of vibrational states, corresponding to different 
excited electronic states, the moving nuclei ATAS spec¬ 
trum of H 2 is very complex, and we will focus on the 
overall structure, interference phenomena in the absorp¬ 
tion lines and LIS. Finally Sec. |VI| concludes the paper. 
Atomic units (h = e = rrif. = Oq = I) are used throughout 
unless indicated otherwise. 


II. THEORY 


The description of a laser pulse propagating though 
a medium in principle requires a full quantum mechani¬ 
cal description of both the field and the medium. In this 
work, however, all pulses considered are far too intense to 
show any non-classical behavior, and the fields are there¬ 
fore treated classically throughout. Even in this semi- 
classical limit the full description of ATAS requires that 
we solve the TDSE and Maxwell’s equations simultane¬ 
ously, which is computationally very demanding |20j . In 
Ref. [in], it is argued that the ATAS intensity spectrum 
of helium found using a single-atom response model, does 
not vary qualitatively from the spectrum found using a 
description including light propagation if the target den¬ 
sity is not too high. A similar conclusion can be drawn 
from the work presented in Ref. [7]. Since our field pa¬ 
rameters are very similar to the ones used in Ref. m, 
and our main focus is a qualitative description of charac¬ 
teristic molecular effects in ATAS, we expect the single¬ 
system response model to be sufficiently accurate for our 
purpose, and we will use it throughout the paper. To cal¬ 
culate how the XUV pulse is modified by propagation in 
ATAS experiments, we follow the approach of Ref. m- 
From Maxwell’s equations we obtain the wave equation 
for the electromagnetic field E(r,t) 
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where c 137 is the speed of light in vacuum and P(r, t) 
is the time-dependent polarization of the target defined 
as the dipole moment density in the target region. We 
consider fields linearly polarized along the z-direction, 
and propagating in the positive a;-direction. We apply 
the single-system response model, thereby reducing the 
interaction between field and matter to a single point in 
space and for the z-components of E(r, t) and P(r,f) we 
obtain 
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The polarization is given in terms of the z-component of 
the single-atom or single-molecule time-dependent dipole 
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moment 


P{t) =n{d{t)), (3) 

where n is the density of atoms or molecules in the sample 
and {d{t)) is the expectation value of the dipole opera¬ 
tor {d{t)) = with the solution of the 

TDSE. The classical field in Eq. ([^ is therefore linked 
to the quantum evolution of the system through P{t). 
We now split the laser field in an incoming field and a 
field generated by the light-matter interaction. We as¬ 
sume that the incoming field propagates unchanged for 
all times 


E{x, t) = Bin (t- + Egen{x, t) . 


(4) 


All modifications of the total field are thus included in 
the generated field Egen{x,t). For the electric field of 
Eq. the wave equation ([^ can be solved analytically. 
The solution for the generated field is (TH] 
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where 0{x) is the heaviside function. From Eq. © we 
see that the generated field has two contributions prop¬ 
agating in opposite directions. 

In the observation point x, the intensity spectrum cor¬ 
responding to the electromagnetic field E{x,t) is ex¬ 
pressed as 
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where Ein{uj) and Egen{uj) are the Fourier components 

Ein{u}) = ^ f ^ j exp{iujt)dt (7) 

1 r°° 

Egen{uj) = — / Egen{x,t)exp{iu}t)dt (8a) 

Ztt J-oo 

/OS 

{d{t)) exp{iujt)dt. (8b) 

oo 


cuin 

c 


Equation (8b I is obtained from Eq. (^8a| by partial in¬ 


tegration, only retaining the solution propagating in the 
positive a;-direction, and finally Eq. (|^ is used. We now 
consider the three terms of Eq. ([^. The first term con¬ 
tains only information about the incoming field, and is 
therefore not of any interest, since we want to investigate 
the modification of the XUV pulse induced by the tar¬ 
get. The two remaining terms contain information about 
this modification. However, since Egen « Fin, the second 
term dominates over the third. Henceforth we therefore 
refer to 


S{u},t) = 2Re[£’*,(a;,r)£;ge„(w,r)], (9) 


as the ATAS spectrum. As seen from Eq. ([^ absorption 
(emission) is described by negative (positive) values of 
S{u},t). In Eq. (|^ we have introduced the variable time 
delay r between the NIR and the XUV pulses in the 
incoming field. It should be noted that the ATAS spec¬ 
trum of Eq. © is not the full intensity spectrum, but 
the modification of the intensity spectrum as compared 
to free propagation of the laser pulses. We refer to the 
full intensity spectrum of ATAS experiments, given by 
Eq. ©, as the ATAS intensity spectrum. For the ATAS 
spectrum to give information about the modification of 
the XUV pulse, it is used that the frequency profile of the 
NIR field is well-separated from the frequency profile of 
the XUV field. Furthermore, the frequency profile of the 
higher-order harmonics of the NIR field must also be well 
separated from the spectral range of the XUV pulse. We 
insure that this is the case by choosing NIR pulses with 
moderate intensities such that only few-photon transi¬ 
tions are induced by these. 

In closing this section we relate our expressions for 
the ATAS spectrum and intensity spectrum to other ex¬ 
pressions used in the literature for similar quantities. In 
Ref. [10] (see also Ref. uni) the response function 


*5(a;)Ref. [ 20 ] = 47rlm[£’*„(a;)d(a;)] 
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is used as a measure of the modification of the ATAS 


intensity spectrum. In Eq. (10), d{uj) is the Fourier com¬ 


ponent d(uj) = ^ /_“ (d(f))e*“*dt. Apart from the sign 
convention, our approach to calculate the ATAS spec¬ 
trum is therefore very similar to the approach of Ref. [20j . 
Another widely used approach to calculate the ATAS in¬ 
tensity spectrum is to use Beers law (see, e.g., Eq. (24) 
in Ref. |2I|) 


|Exuv(E,a;)p =|£'xuv(0, w)| 
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where Exxsy{L,uj) is the Fourier component of the XUV 
field when it has propagated the distance L into a homo¬ 
geneous target medium, and P{oj) is the Fourier compo¬ 
nent of the polarization. For weak absorption it is accu¬ 
rate to expand Eq. 0 to first-order in the argument of 
the exponential function, and we find that 

|Exuv(A,w)p |Uxuv(0,w)P + LS{uj,t). (12) 


This means that both here, and in the theory of Ref. m, 
the ATAS spectrum S{u}, r) of Eq. © is a good measure 
of the modification of the XUV intensity spectrum when 
the generated field is weak {Egen « Ejn). Some times 
the cross section rather than the intensity modification 
is reported in the literature. This is, e.g., the case in 
Refs. |7l[22|. In Ref. |20| it is shown that the absorption 
cross section cr(w, r) is directly related to S{uj, t) through 
cr(w,r) = -2/(n\Ein(uj)\^)S(uj,T). 
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III. RESULTS FOR ID H2+ BASED ON TDSE 
CALCULATIONS 


In this section we present the ATAS spectru m of 
Eq. ([^ for two different ID models of H 2 ^. In Sec. Ill A 
we present the ATAS spectrum calculated in the hxed 
nuclei model of H 2 . For this system we recognize many 
features from atomic spectra previously reported in the 


literature In Sec. |IIIB| we present results 

for the moving nuclei model of 9'iid we find that 
dissociation processes completely change the spectrum. 

We define the z-component of the electric field E{t) = 
-dtA{t) through the ID vector potential 


A{t) = Aq exp 


{t-tcf 


T‘^14 


cos [uj{t - tc)] , 


27r 


T=N,Tc = N, — , (13) 

OJ 


(TOp + l)/mp, rup = 1.836x10^ is the proton mass, and the 
softning parameter oq = 1.508 for the Coulomb singular¬ 
ity is chosen to produce the exact three-dimensional BO 
electronic energy at i?o- A box size of \z\ < 100, Az = 0.39 
and time steps of At = 5x 10“^ ensured converged results. 
At the box boundaries a complex absorbing potential was 
used to avoid unphysical reflections (see Ref. [23] for de¬ 
tails) . 

Following excitation by the XUV pulse, an oscillating 
dipole moment is established in the target molecule. In 
experiments the time-dependent dipole moment dephases 
during a time Tq due to collisional broadening and fi¬ 
nite detector resolution m- To include these effects we 
use a window function W{t) to damp (d{t)) on a time 
scale, that produces spectra comparable to experimen¬ 
tal results. Further, we ensure that the dephasing time 
is long enough for the qualitative features of the ATAS 
spectrum not to vary with Tq. For a window function 


where Aq = EqIuj. Eq relates to the intensity as I = 
l-Fop, with one atomic unit of intensity equal to 3.51 x 
IQi® W/cm^. In Eq. (13) tc is the center of the pulse 
and Lo is the carrier frequency. The full width at half 
maximum (FWHM) duration Tfwhm is related to T used 
in Eq. (13) by Tfwhm = ’yin(2)r. For the calculations 
reported in this section, we use a 700 nm NIR pulse with 
T = 14.0 fs [Ac = 6 in Eq. ( [I^ ] and a 50 nm XUV pulse 
with T = 330 as [Ac = 2]. The intensities of the pulses 
are /xuv = 5 x 10^ W/cm^ and Inir = 2 x 10^^ W/cm^, 
respectively. At these intensities the NIR pulse alone 
cannot excite the molecule from its ground state to the 
spectral region of the XUV pulse, as many-photon (> 3) 
processes are greatly suppressed. To propagate solutions 
of the TDSE for the two systems we use a fast Fourier 
transform (FFT) split-step operator algorithm [53], as 
described in Refs. 1531 [55]. 


A. Fixed nuclei model 


We consider a model for H 2 ^ of reduced dimensionality 
with the fixed internuclear distance Rq = 2.07 which is 
the expectation value of the internuclear distance in the 
ground state. The TDSE in the dipole approximation 
and velocity gauge reads 

RqA) = A(t)T(z; Rq, t) (14) 

with the Hamiltonian 

H{t) = Te + UeN(^; i?o) + V^{z-Rq) + Vl{t) , (15) 


where z is the electronic coordinate measured with re¬ 
spect to the center-of-mass of the nuclei. The com¬ 
ponents of the Hamiltonian in Eq. (15) are Te = 
-(l/2/i)5Vaz2, UeN(Ai?) = -l/y(^A-R/2)2T^ _ 

l/V(2+i?/2)2 + ao, V^{R) = IjR, Vi{t) = -ipA{t)djdz, 
= 2mp/(2mp + 1) is the reduced electron mass, /? = 


W{t) = 


i{l + cos[^]} 


(t < r) 

(r < t < T + To) 

(t > r + To), 


(16) 


where t is measured with respect to the center of the 
NIR pulse, we find that Tq = 3Tnir - 42 fs is suitable. In 
Fig. we present the ATAS spectrum [Eq. (§] 

calculated for the fixed nuclei model of H 2 ^, with delays 
of the XUV pulse with respect to the NIR pulse varying 
from -28 fs to 14 fs. For large positive delays, the NIR 
pulse arrives at the target well before the XUV pulse. 
Because the relatively weak NIR pulse leaves the system 
in its held-free ground state, it will not have any effect on 
the ATAS spectrum, when it arrives well before the XUV 
pulse. The spectrum therefore shows clean absorption 
lines without side-bands or other modifications for delays 
larger than ~ 8 fs. The position of the four absorption 
lines in Fig.[2at approximately 10 eV, 21.5 eV, 25.5 eV, 
and 27 eV, correspond to resonances between the field- 
free ground state and states of opposite symmetry under 
the reflection z -*■ -z with energies 9.94 eV, 21.48 eV, 
25.29 eV and 26.94 eV above the ground state energy, 
respectively. 

For large negative delays, the XUV pulse arrives well 
before the NIR pulse. The NIR pulse therefore arrives 
when an oscillating dipole moment has already been in¬ 
duced. Even though the NIR pulse does not excite ad¬ 
ditional population from the ground state to the states 
reached by the XUV field, the NIR pulse will modify 
(d{t)). One effect of this modification is seen around 
the absorption peak at -- 21.5 eV in Fig. [^ For delays 

more negative than-5 fs, we see that the absorption 

line has several sidebands. Such sidebands have previ¬ 
ously been observed in atomic systems (see, e.g.. Ref. m 
and references therein). Another effect introduced by the 
NIR field is a periodic modulation of the absorption line 
at ~ 25.5 eV in Fig. (Enlarged in the upper panel of 
Fig. 1^. This effect extends from large negative delays 
and all the way to the region where the two pulses over- 
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XUV first XIR first 



Delay (fs) 

Figure 1. (Color online) ATAS spectrum [Eq. ([^] 

calculated for the fixed nuclei model of H 2 ^, as a function of 
the photon energy, for a range of delays between the XUV 
and NIR pulses. The ATAS spectrum is normalized such that 
maximum absorption is at S(ai,T) = -1.0. For photon en¬ 
ergies between 14.0 eV and 20.3 eV the ATAS spectrum has 
been multiplied by a constant factor of 50 to highlight the LIS 
doublet in that region (see text). The white line in the upper 
right part of the figure indicates the intensity profile of the 
XUV pulse. Pulse parameters: Txuv = 330 as, Tnir = 14.0 
fs, Axuv = 50 nm, Anir = 700 nm, 7xuv = 5 x 10^ W/cm^, 
/nir = 2 X 10^^ W/cm^ 


lap in time, and is a consequence of quantum interference 
between populated excited states of the same symmetry 
m [12]. For large negative delays the XUV and NIR 
pulses are separated in time, and the role of the XUV 
pulse is to populate the excited states at time r while 
the NIR pulse modifies the populations and coherences 
at a later time tg. In Appendix [A] we describe a (three- 
level) model which explains and reproduces the periodic 
modulation of the absorption line at ~ 25.5 eV. A result 
of this treatment is that oscillations of this type always 
have a period of ~ T^^^/2 with the period of a NIR 
pulse cycle. 

The region -I5fs < r < 5fs, where the two pulses over¬ 


lap in time, is dominated by intense coupling between 
states. As a consequence, absorption, and to a smaller 
extent emission, appears at certain frequencies in this re¬ 
gion only. Here we will focus on two features both visible 
in Fig. 1^ The Autler-Townes splitting of the ab¬ 
sorption line at ~ 21.5 eV, and the LIS [10] doublet with 
spectral features at ~ 19.4 eV and ~ 15.9 eV correspond¬ 
ing to energies one photon energy above and one photon 
energy below the even parity state at 17.65 eV. When 
the two pulses overlap in time, we see a clear broadening 
and a splitting of the absorption line at ~ 21.5 eV. Such 
features originate from coupling between bound states 
and have been observed both experimentally m and in 
calculations m- As discussed in Refs. mm, the nature 
of the Autler-Townes splitting can be understood in few- 
state models. LIS are signatures of population transfer, 
by two-photon processes, from the ground state to an 
excited electronic state of the same symmetry, under the 
reflection z ^ -z. In spectroscopy such a state is often 
referred to as a dark state. Two different processes, one 
involving absorption of a NIR photon and one involving 
emission of a NIR photon, can populate the same dark 
state. As a result LIS always appear in pairs with an 
energy separation of 2a;NiR between the individual ab¬ 
sorption features. The two intermediate states absorbing 
XUV light are always in two NIR-photon resonance re¬ 
sulting in strong coupling between the states m- This 
coupling leads to a rj^^^/2 periodic modulation of the 
LIS. In Ref. [10] it has been shown that specific LIS in he¬ 
lium can be reproduced using only three states in the cal¬ 
culation of the ATAS spectrum. For H 2 ^ with fixed nuclei 
we set up a similar three-level model including the ground 
state, a dark state 17.65 eV above the ground state, and a 
dipole allowed state 21.48 eV above the ground state. We 
will refer to these field-free states as (f)g{z;Ro) (ground 
state), (dark state), and cj)f,{z-,Ro) (dipole al¬ 

lowed, excited state) and their corresponding energies as 
Eg = 0, Ed, and E^- The wave function of a general state 
in the three-level model is given by 

4'(z;i?oR) = Cg(t)(j)g(z-Ro) + Cd(t)(j)d(z; 

+ (17) 


We find that the three-level model reproduces the LIS 
shown in Fig. at Ud±WNiR well. For the state (17) the 
expectation value of the dipole operator is well approxi¬ 
mated by 


(d(t))^2Re[dgeCe(t)e-*^'(‘-")], (18) 

since the condition Cg{t) 1 » |ce(t)| > \cd{t)\ holds for 
field parameters used in this work. 

In the pursuit of a simple analytical theory, capable 
of explaining the characteristics of LIS, we now simplify 
our three-level model further. This analytical result will 
be needed for understanding the moving nuclei results 
of Sec. jVj The spectral widths of LIS are determined by 
the duration of the NIR pulse. This indicates that the 
temporal region where the XUV field is on only has a 
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small direct effect on the ATAS spectrum (The Fourier 
transform of the dipole moment in this region is van¬ 
ishing when compared to the contribution from the re¬ 
maining temporal region of the NIR field). The main 
contribution from the temporal region where the XUV 
and NIR pulses are overlapping is therefore to populate 
the states (j)d{z-,Ra) (two-photon process) and (j)eiz-,Ro) 
(one-photon process). When the shorter XUV pulse has 
died out, the population transfer from (j)g{z;Ro) to the 
excited states will stop, and the dynamics will be dom¬ 
inated by coupling between (pdiz^Ro) and (j)e{z]RQ) in¬ 
duced by the NIR field. When we consider instants of 
time after the XUV pulse {t > r) but before the end of 
the NIR pulse, the dynamics of the system is approxi¬ 
mately governed by the equations 


Cg{t) = 0 (19) 

Cd{t) = (20) 

Ceit) = (21) 


In Eqs. (|19[)-(21 ) the multi-photon NIR-coupling between 
the ground state and excited states has been neglected 
since such multi-photon processes are greatly suppressed 
for the held parameters used in this work. If we for 
simplicity assume that the NIR held is given by the 
monochromatic held E{t) = U'^jj^sin(wNiRt), we have 


LIS established in this section will be the starting point 
of an investigation of LIS in molecular systems in Sec.[V| 
Figure and the associated discussion recapitulate 
many of the physical processes already investigated in 
atoms. The detailed spectral and temporal information 
about electron dynamics extracted from Fig.[l]is a conse¬ 
quence of the long dipole response of the system, allow¬ 
ing for spectral resolution much better than the Fourier 
limit of the XUV pulse. With the phenomena discussed 
in this section in mind, we are ready to explore the effect 
of molecular vibrations. 


B. Moving nuclei model 

In this section we present the ATAS spectrum of 
Eq.(l|) for the moving nuclei model of H 2 ^ using the same 
held parameters as in Sec. |III A] We consider a simplihed 
model for H 2 ^ with reduced dimensionality that includes 
only the dimension that is aligned with the linearly po¬ 
larized laser pulse [HHII]- Within this model, electronic 
and nuclear degrees of freedom are treated exactly. The 
center-of-mass motion of the molecule can be separated, 
such that the TDSE for the relative motion in the dipole 
approximation and velocity gauge reads 

idt'i>{z,R,t) = H{t)'i>{z,R,t) (24) 


Vdeit) = V,d{t) = - e-™*]. (22) 

For the held parameters used in Fig. [^the detuning be¬ 
tween (l)d{z', Ro) and (j)e{z] Rq) is large {Eg-Ed = 3.83 eV, 
i^NiR = 1-77 eV) and the population transfer is therefore 
low. An approximate analytic expression for Ce{t) can 
therefore be found using perturbation theory. First or¬ 
der perturbation theory can reproduce the characteristics 
of the LIS and is desirable since the result is relatively 
simple 




TT’O J 

^NIR^rfe 

2i 


^i{Ec-Ea+LUNiii)x{t-T) _ ^ 


Ee- Ed + WnIR 


^i(E^-Ed-LOKm)x(t-T) _ ^ 
-f 

Ee - Ed- wnir 


Cd{T). (23) 


If we insert Eq. (23) into Eq. 


we see that {d{t)) 


will contain terms oscillating with frequencies Ed iWNiR. 
S{oj) [see Eq. (|^] found using this dipole moment, will 
therefore contain spectral features at Ed ± wnir - the LIS 
seen in Fig. [T] It should be noted that the part of the 
time-dependent dipole moment responsible for the LIS 
depends critically on the initial population of the dark 
state, which is in compliance with the two-photon na¬ 
ture of LIS. From Eqs. (23) and (18) we conclude that 
the NIR-coupling of the excited states in our three-level 
model gives rise to LIS. The simple analytical model for 


with the Hamiltonian 


H{t) = Te + Tn + Ve^{z,R) + V^{R) + Ui(t), (25) 

where the wave function now depends on the internu- 
clear distance R and the electronic coordinate z mea¬ 
sured with respect to the center-of-mass of the nu¬ 
clei. The components of the Hamiltonian in Eq. (25) 


were defined after Eq. (15), except the nuclear ki¬ 
netic energy Tjv = IdR? and I 4 n(- 2 ,7?) = 

-l/y/(z-i?/2)2 + a{R) - l/^/(z + i?/2)2 + a(i?), where 
the softning parameter a(i?) for the Coulomb singular¬ 
ity is chosen to produce the exact three-dimensional IsUg 
BO potential energy curve [231 HH]- A box of |z| < 100, 
Az = 0.39 in the electronic coordinate together with a 
box of nuclear coordinates of 0 < i? < 20, Ai? = 7.8 x 10“^ 
and time steps of At = 5x 10~^ ensured converged results. 
At the box boundaries complex absorbing potentials were 
used in both electronic and nuclear coordinates to avoid 
unphysical reflections (see Refs. [231|2S] for details). 

The ATAS spectrum is shown in Fig. It is clear from 
Fig. I^that the relatively sharp lines we saw for the hxed 
nuclei model [Fig. have disappeared. Two very broad 
absorption lines corresponding to population transfer be¬ 
tween the ground and the first excited state (absorption 
feature at ~ II eV) and between the ground state and the 
third and fifth excited states (merged absorption feature 
at ~ 23 eV) dominate the spectrum. The broad features 
in Fig. [^persist as r ^ ±oo, and are therefore not a result 
of the NIR pulse. The only effect of the NIR pulse in Fig. 
[^is the small oscillations with a period of T^®/2 clearly 
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Figure 2. (Color online) Normalized ATAS spectrum S(aj, r) 
[Eq. calculated using the moving nuclei model [Eqs.([M| 
and (25l], for ID H 2 ^, as a function of the photon energy, 
for a range of delays between the XUV and NIR pulses. The 
white line to the right shows the intensity prohle of the XUV 
pulse. The pulse parameters are as in Fig. 


present between ~ 14-28 eV for delays -5 fs < r < 5 fs, 
corresponding to the region where the pulses overlap the 
most. In Sec. |1V| we investigate the molecular dynamics 
resp onsibl e for the broad absorption lines in Fig. In 
Sec. |lV B we also briefly discuss the origin of the T^^/2 
oscillations in Fig. 


IV. ANALYSIS AND DISCUSSION OF H2+ 
RESULTS 

In this section we set up a simplified model to explain 
the broad absorption features in the ATAS spectrum of 
the moving-nuclei model of H 2 ^. To isolate the effect we 
first neglect the NIR pulse, whose impact we found to be 
small [Fig. [^. The analysis will show that the dissocia¬ 
tive nature of the excited electronic states in H 2 ^ is re¬ 
sponsible for the broad absorption features in the ATAS 
spectrum and will provide a link between the spectral 


width of the ATAS spectrum and the time it takes for 
the part of the nuclear wave packet found in an excited 
electronic state to dissociate. 


A. A-surface model 

In the fV-surface model the wave function of the 
molecule is expanded in the lowest N electronic BO 
states, (j)j{z]R)-. 


N 

= (26) 
j=i 

where Gj{R,t) are the corresponding nuclear wave pack¬ 
ets and z denotes the electronic coordinates. In the N- 


surface model of Eq. (26), the time-dependent expecta¬ 
tion value of the z-component of the dipole operator d, 
entering the evaluation of the ATAS spectrum as dis¬ 
cussed in Sec. |nj is given by 

{d{t)) = / dR J dz(z, R,t)d'^{z, R,t) 

N f. 

= J dRGKR,t)G,{RR)d:y{R), (27) 


where 

d:y{R) = f dz^:iz-,R)dMz;R) (28) 

is the electronic dipole moment function. The evolution 
of the nuclear wave packet Gi{R, t) on the energy-surface 
Ee\,i{R), corresponding to the Fth electronic state, is 
governed by the equation 


N 

idtG,{R,t) = +Uy)G,(R,t), 


(29) 


obtained by inserting the iV-surface expansion (26) in the 
TDSE [Eq. (24)] and neglecting terms containing deriva¬ 
tives of the electronic states (j)j{z,R) with respect to R. 
In Eq. (|2|) = -^^ + Eei,,(R) where E.i^R) 

is the BO potential energy curve obtained by solving 
the electronic problem at fixed R. Due to the trunca¬ 
tion in the number of electronic states, the A-surface 
expansion (29) is not gauge invariant, and the dynamics 
are only correctly described in the length gauge [29] . In 
Eq. (29) we therefore have that Vij = -fdLcdfj, where 


Plg = {Zi + Z 2 + Ml + M 2 )/Mtot, with Zi and Z 2 the nu¬ 
clear charges. Mi and M 2 the nuclear masses and Mtot 
the total mass of the molecule. Note that an atomic ver¬ 
sion of the A-surface model was used for He in Ref. [12j . 

For a further analysis of the broad absorption features 
in Fig.[^ we expand the nuclear wave packets Gi{R,t) in 
the vibrational eigenstates Xi,k{R) corresponding to the 
j’th electronic state 


G,{R,t) = ^dE,,k c,,fc(t)x,,fc(i?)e-*^-'‘*, (30) 
























Nc = 3 


where Ei^k are the eigenenergies of Xi,k{R)- The set of 
Xi,k{Rys for a specific i can be continuum states, or a 
mixture of bound and continuum states. The continuum 
states are energy normalized, such that 

f dRxlp{R)x^AR) = - E^,k)■ (31) 


Using the expansion (30), the time-dependent dipole mo¬ 


ment of Eq. (27) reads 

= E ^ dEj^p c*At)cjAt) {x^,k\dfj IxjA 

AEi,k-Ej^p)t ^ ( 32 ) 

where 

{x^Mdt \Xj,p) = f dR xlk{R)Xj,p{R)dfAE)- (33) 

Often the function dfj{R) varies slowly with R, as long 
as R is varied in a region where G*{R,t)Gj{R,t) is non¬ 
zero. In such cases it is a good approximation to let 
dij{R) = dijiRo) and we obtain 

N 


{d-{t)) = E dEi^k ^ dEj^p c*At)cjAt)dfj{Ro) 


(Xi,k\X3.p)e 




(34) 


where the subscript indicates that the variation in R is 
neglected for the function In Eq. (34) {xi,k\Xj,p) 

is the Franck-Condon overlap between the states XixiR) 
and XjAR)- 

The intensity of the XUV pulse is generally very low 
in ATAS, and it is reasonable to assume that the ground 
state population remains unchanged during the pulse, 
i.e., ci_o(l) - Ij and that the expansion coefficients Ci^k{t) 
for vibrational states, corresponding to the BO surfaces 
of excited electronic states, can be found using first-order 
perturbation theory 

(^) = f dtEinit) (xi,k\df^i |xi,o) 

(35) 


B. A-surface model for H 2 ^ 

In H 2 ^ all excited BO surfaces are dissociative, and 
we therefore expect {d{t)) to be damped when the ex¬ 
cited nuclear wave packets dissociate [see Eq. (27)] This 
damping leads to a broadening of the ATAS spectrum. 
For H 2 we found that the moving nuclei ATAS spectrum 
of Fig. 1^ can be reproduced to a very high accuracy using 
an A-surface expansion [Eq. (26)] of the total wave func¬ 
tion followed by propagation of the nuclear wave func¬ 
tions Gi{R,t) in accordance with Eq. (29). The nuclear 
wave packets were propagated using a split step algo¬ 
rithm, A = 6 in Eq. (29), 0 < i? < 20, Ai? = 7.8 x 10"^ 



10 20 30 40 10 20 30 40 

Photon energy (eV) 

Figure 3. (Color online) ATAS spectrum [Eq. §] 

(solid, black lines) and field intensity (dashed, red lines) as 
functions of the photon energy. In the left column pulses 
with Ac = 1 in Eq. (131 have been used in the calculations, 
while pulses with Ac = 3 have been used in the right column. 
The central carrier frequencies of the pulses used in the cal¬ 
culations are 6 eV for (a) and (b), 10 eV for (c) and (d), 20 
eV for (e) and (f), 30 eV for (g) and (h), and 50 eV for (i) 
and (j). 


and time steps of At = 5x 10~^ ensured converged results. 
The fact that ATAS spectra can be precisely determined 
in the A-surface model is of importance for studies of 
ATAS in larger systems which can not be treated by 
direct solution of the TDSE. In Sec. |V] we treat H 2 in 
three dimensions using the A-surface model. A very use¬ 
ful property of the A-surface expansion is that we can 
subtract specific surfaces from the expansion during the 
propagation, and in this way determine their influence 
on the ATAS spectrum. In Fig. we saw oscillations 
with a period of T^^^/2 in the region 14-28 eV. By per¬ 
forming separate calculations, we have found that these 
oscillations are due to two-photon processes involving the 
electronic even-parity states (p 3 {z\R) and (p 3 {z;R). 

In H 2 ^ the coupling from the ground state to the 
first excited electronic state is several orders of magni¬ 
tude larger than the coupling to higher excited electronic 
states when evaluated near R = Rq. It is therefore ex¬ 
pected that pulses with broad frequency distributions pri¬ 
marily populate the first excited state. This leads to an 
interesting effect for ultrashort pulses; the shape of the 
ATAS spectra shown in Fig. for pulses with Ac = I in 
Eq. (13) and carrier frequencies corresponding to energies 


higher than 20 eV, turn out to be remarkably unchanged 
as the carrier frequency is varied. As can be seen from 
Fig-i this is not the case for longer pulses. The FWHM 
of the three pulses with Ac = 1 and central carrier fre¬ 
quencies corresponding to 20 eV, 30 eV and 50 eV are 
2.67 eV, 2.67eV and 2.64 eV, respectively. This indi- 
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cates that the ATAS spectra for these ultrashort pulses 
are only determined by the system, and that only the 
ground state (j)i{z]R) and the lowest excited electronic 
state (j) 2 {z]R) are substantially populated. In the case 
of an ultrashort pulse the dipole signal is therefore ex¬ 
pected to be well approximated by a simplified version of 
Eq. PI 


(di^^(t))=2Re| -z J dE2,kU{E2,k)df^dRo) 


l(X2,fc|xi.o)l e 




(36) 


UiE2,k) = f 

where the superscript ’(!)’ indicates that we have used 
first-order perturbation theory [Eq. (35)1 to determine 
the coefficients in Eq. ( p| . Equation (|36| is a good ap¬ 
proximation for the time-dependent dipole moment as 
confirmed by the results in Fig.[^ which will be discussed 
later. 

After the pulse is over U{E 2 ,k) does not vary with 
time. In H 2 ^ the FWHM AE of |(x 2 ,fc|xi,o)l^ is ~ 2.5 eV. 
For pulses of duration much less than = 1.7 fs we can 
therefore neglect the variation of U{E 2 ^k) with respect to 
E^.k] U{E 2 ,k) = U. For such short pulses we have 



Time (fs) 


Figure 4. (Color online) Time-dependent dipole moment 
from the full TDSE calculation (solid, black line) and the 
approximation of Eq. (361 (dashed, red line) calculated 
for the ID model of H 2 ^ and a single-cycle 30-nm pulse. The 
intensity of the pulse is 10^ W/cm^. The light blue dotted 
line shows a Gaussian fit to the norm of the extrema in the 
damped dipole moment. The hrst data point included in the 
fit is the extremum at t 0.7 fs. The XUV pulse is centered 
at t = 0.5 fs. 


(di^^(t)) = 2Re 


-iUdf^,{Ro)E^^’°* f dE2,k 

{t» E) (37) 


|(X2,fe|xi,o)|^-*''^-''‘ 


the 


pulse. 

inverse Fourier 
2 


where tc is the center of 

/c^F 2 ,fc|(x 2 .fe|xi.o)f is the 

transform of the Frank-Condon overlaps |(x2,fc|xi,o)| 
with respect to the Fourier variable E 2 .k- It turns 
out that |(x2,/c|xi,o)|^ is well approximated by a 
Gaussian function D{E 2 ,k - E 2 {Ro)) centered around 
E 2 .k - E 2 {Ro) = 0, where E 2 {Ro) is the value of the first 
excited BO surface, evaluated at Rq = 2.07. Substitution 


with the variable E 2 .k = E 2 .k - E 2 {Ro) in Eq.(37) gives 


(d[^\t)) = 2Re|^ - 

J dE2^kD{E2,k) 


-iE2,k't 


{t » E) (38) 


83 as, and we therefore expect the corresponding dipole 
signal to be well described by Eq. (38), when the pulse is 
over. A Gaussian fit to the absolute values of the peaks of 
the dipole signal, in the region where the pulse has died 
out (more than 165 as, ie., 2 FWHM durations, after the 
central time of the XUV attosecond pulse) is shown in 
Fig.H and we see that the damping of the dipole signal is 
well described by a Gaussian function. We also see that 
the dipole signal is oscillating with a period close to the 
expected period of 420 as. 

Assuming that the Fourier transform of the incoming 
field is constant in the region where the ATAS spec¬ 
trum is non-zero (see Fig. [^, the ATAS spectrum 
of Eq. ([^ is proportional to the Fourier transform of 
the time-dependent dipole moment [neglecting the lin¬ 
ear frequency dependence from Eq. (8bl]. For ultra- 
short pulses the ATAS spectrum of H 2 '** is therefore well- 
approximated by a Gaussian function, and the half-life 
ti /2 of the envelope of {d{t)) is related to the FWHM 
wfwhm of the ATAS spectrum through the relation 


Since the Fourier transform of a Gaussian is also a 
Gaussian, we can now conclude that the time-dependent 
dipole moment of (38), after the pulse is over, is oscil¬ 
lating with frequency Ei q - E 2 {Ro) and is damped by a 
Gaussian function. The period of oscillation correspond¬ 
ing to the energy Ei^o-E 2 {Ro) is ~ 420 as. Figure ^shows 
the time-dependent dipole moment from H 2 ^ calculated 
using a 30 nm pulse with an intensity of 10^ W/cm^ and 
Nc= 1 in Eq. (13). The FWHM duration of this pulse is 


^ 1/2 - 


41n(2) 

WfWHM 


(39) 


From the three ultrashort pulses discussed earlier (left 
panel in Fig. with wfwhm - 2.65 eV we find that 
tii 2 - 690 as, which is indicated on Fig. The time 
ti /2 indicates the duration it takes the time-dependent 
dipole-moment to decay and it therefore sets the time 
scale in which the nuclear wave packet on the excited 
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state G 2 {R,t) travels to a region where the overlap with 
the ground-state wave packet Gi{R,t) is significantly de¬ 
creased. The value of ti /2 shows that this dynamics takes 
place on a subfemtosecond time scale. This very fast 
nuclear motion explains the extremely broad absorption 
features in the spectrum shown in Fig. 


V. Ha IN THREE DIMENSIONS 


In all BO surfaces corresponding to excited elec¬ 
tronic states lead to dissociation. In molecules in general 
there are also excited BO surfaces supporting bound nu¬ 
clear motion. To have a more general description of nu¬ 
clear motion effects in ATAS, we therefore consider the 
H 2 molecule which is sufficiently simple to allow for an 
accurate numerical treatment, but still show the charac¬ 
teristics of excited BO surfaces with bound nuclear mo¬ 
tion. 

Inspired by the good performance of the A^-surface 
model in ID H 2 ^, we represent the total wave function of 
3D H 2 in the A^-surface expansion [Eq. (26)]. The nuclear 
wave functions Gi{R,t) are propagated in accordance 
with Eq. (29) using a EET split-step algorithm. BO sur¬ 


faces and dipole transition moments used in the propa¬ 
gation are given in the literature m- We must ensure 
that the basis used in the expansion of Eq. (26) is large 
enough to contain all physics relevant to the ATAS spec¬ 
trum. To reduce the number of basis states N needed, 
we choose the intensity of the dressing field low enough 
that processes involving more than two low-energy pho¬ 
tons are suppressed. Calculations of for H 2 with 

fixed nuclei did not show any signs of 3-NIR-photon or 4- 
NIR-photon processes for the field parameters we use in 
this section. To further reduce N, we use a 1600 nm IR 
dressing pulse. Under these criteria it is sufficient to use 
Af = 10 BO surfaces in the expansion. The propagation 
of Gi{R,t) with 0 < i? < 20, AR = 7.8 x 10”^ and time 
steps of At = 5 X 10“^ ensured convergenced results. To 
highlight the effect of nuclear dynamics on ATAS spectra 
we also preform a fixed nuclei calculation on 3D H 2 using 
the same number of electronic states as described above. 
In the fixed nuclei calculation we use an internuclear sep¬ 
aration of i? = 1.393, corresponding to the separation at 
minimum energy in the ground state BO curve. 

In Eig.j^we present the ATAS spectrum calcu¬ 

lated for H 2 with delays of the XUV pulse with respect to 
the IR pulse varying from -84 fs to 10 fs. Eigure Fig.j^a) 
shows the fixed nuclei results, while Fig. ^b) shows the 
ATAS spectrum, when the nuclei are allowed to move. 
The field parameters are the same in the two calcula¬ 
tions. The IR pulse is a 1600 nm pulse with a FWHM of 
8.83 fs [Nc = 2 in Eq. and the XUV pulse is an 80 
nm pulse with a EWHM of 440 as {Nc = 2). The intensi¬ 
ties of the pulses are Ar = 10^^ W/cm^ and Ixuv = 5x 10^ 
W/cm^, respectively. The 11th (singlet) BO surface in 
H 2 has an energy minimum of ^ 14.5 eV above the ground 
state energy EH, and the photon energy of the IR field is 
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Figure 5. (Color online) Normalized ATAS spectrum ^(a;, r) 
[Eq. ([^] calculated for 3D H 2 . (a) for fixed nuclei, (b) for 
moving nuclei. The shaded areas indicate photon energy re¬ 
gions where the spectra are expected to be influenced by the 
finite electronic basis (see text). The white line to the right 
indicates the intensity profile of the XUV pulse. The dots 
and triangles to the very right in (b) are located at the field- 
free energies of the vibrational states corresponding to the 
BO curves B“ (black dots) and C“ (black triangles). Pulse 
parameters: Txuv = 530 as, Tir = 10.7 fs, Axuv = 80 nm. 
Air = 1600 nm, 7xuv = 5 x 10^ W/cm^, Ar = lO'^^ W/cm^. 


0.775 eV. In Fig. 5 regions of photon energies above 14.5 
eV -2 X 0.775 eV = 12.95 eV are shaded indicating that 
we expect structures in this region of the ATAS spectrum 
not to be well described in the basis of 10 BO surfaces. In 
H 2 we, in addition to electronic dynamics, want to follow 
nuclear dynamics occurring on a much longer timescale. 
We therefore use a longer dephasing time for the time- 
dependent dipole moment than for the case of H 2 with 
fixed nuclei. For the calculations in this section we use 
the window function (16) with Tq = 48 fs, but shift the 
window function such that {d{t)) is undamped the first 
73 fs after the end of the XUV pulse. 


From Fig. [^we immediately see that the ATAS spec¬ 
trum of H 2 is highly affected by the nuclear dynamics of 
the system. The effect of nuclear motion on the ATAS 
spectrum in H 2 is, however, very different from what 
we saw in H 2 . In the fixed nuclei ATAS spectrum of 
Fig. I^a) we only see a single absorption line in the un¬ 
shaded region. For H 2 with a fixed internuclear sepa¬ 
ration of i?o = 1.393, the excited states considered have 
energies higher than the lowest vibrational state in the 
corresponding BO curve for the moving nuclei system. It 
is therefore expected, that we can go to higher energies in 

























11 


the ATAS spectrum without expanding the basis of the 
calculation. For energies higher than 12.95 eV additional 
absorption lines are visible in the fixed nuclei ATAS spec¬ 
trum [one is included in the shaded area of Fig. [^b)]. 
These absorption lines and the absorption line at 12.8 
eV all behave very much like the absorption lines in the 
fixed-nuclei spectrum of H 2 ^; we observe interference of 
the type discussed in Appendix A between states in two 
NIR photon resonance, Autler- Townes splittings of ab¬ 
sorption lines etc. 

The moving nuclei ATAS spectrum of H 2 shown in 
Fig. j^b) is, on the other hand, very different from both 
the fixed nuclei ATAS spectrum of H 2 , and from the 
moving nuclei ATAS spectrum of H 2 ^ (Fig. [^. Even 
though only three singlet BO curves (B“, C“ and EF®) 
are within the energy range considered in Fig. [^b), we 
see rich structures in the ATAS spectrum. For the pos¬ 
itive delays to the outermost right in Fig. ib) we see 
that these structures converge towards absorption lines 
corresponding to population of the field-free vibrational 
states of the B“ and C“ surfaces. When the delay be¬ 
tween the XUV pulse and the NIR pulse is varied, we 
see that oscillations occur in the spectrum on two dif¬ 
ferent time scales. Throughout the spectrum there are 
oscillations on the time scale of half an IR period (~ 2.7 
fs), and in the absorption lines corresponding to vibra¬ 
tional states of the B“ surface [indicated by black dots in 
Fig-[5];b)], we see additional oscillations on a much longer 
time scale. Introducing the Fourier transform of 
[Eq. (§] 

1 r°° ~ 

.F[5'(a;,T)] = — / exp{iujrT)dT, (40) 

Ztt J—00 


we are able to analyze the oscillations of Fig. [^b) closer. 
The Fourier transform Eq. (40) of the delay-dependent 
ATAS spectrum is shown in Fig. I^a) for low and in 
Fig. I^b) for higher frequencies Wr- The low frequency 
structure in the Fourier transform of the ATAS spectrum 
is a consequence of nuclear dynamics in the molecule. In 
the Wsurface expansion (26), we can think of the XUV 
pulse as generating nuclear wave packets on BO surfaces 
corresponding to different excited electronic states. If the 
BO surfaces were harmonic, the wave packets would os¬ 
cillate periodically with a period of 2t: jujh, where Wh is 
the frequency of the harmonic oscillator. However, the 
BO curves of H 2 are not perfectly harmonic, and we must 
account for their anharmonicity. We do this by defining 
a local vibrational period for the n’th vibrational state 
as 


J.I 0 C ^ 1 / 27r ^ 27r \ 

2 y ^^n,n+l / 


(n>l), (41) 



Figur e 6. (Color online) Absolute value of the Fourier trans¬ 
form (401. (a) for relatively low frequencies and (b) for higher 
frequencies. The dashed line in (a) indicates the local energy 
separation between vibrational states in the B“ curve cor¬ 
responding to the local vibrational period of Eq. ( |41[ ). The 
lines in (b) are Ei - uir [see Eqs. (40) and (42)] when Ei is 
either the energy of one of the seven lowest vibrational states 
of the C“ surface (dashed lines), one of the five lowest vibra¬ 
tional states of the B’“ surface (dash-dotted lines), or one of 
the three lowest vibrational states of the D“ surface (dotted 
lines). See Ref. [30] for information on the energy levels. 


the energy separation of the vibrational states in the C“ 
curve around 0.3 eV. 

In Fig. I^b) we present the higher frequency compo¬ 
nents of the Fourier transform of the delay-dependent 
ATAS spectrum. Note that the structures of Fig. |^b) 
are centered around the energy 1.55 eV of two 1600nm 
IR photons. Much of the structure in Fig. I^b) can be 
explained as interference between different paths to the 
same final state, just as the interference in the absorp¬ 
tion lines for fixed-nuclei H 2 ’^ (see Appendix[A|). For both 
fixed and moving nuclei the oscillation has a period of 
approximately half a NIR or IR field cycle. We now 
focus on one absorption peak in Fig. and denote its 
energy by w. If the field-free state corresponding to this 
absorption peak couples to a state with energy Ei, the 
absorption peak at energy w will contain a term propor¬ 
tional to cos[(a;-Ei)r], neglecting a phase. We now take 
the Fourier transformation defined in Eq. (40), but with 
S{uj, t) cos[(w - Ei)T], and see that we only obtain a 
nonvanishing contribution when 


Ei ±uJt = uj. 


(42) 


where AuJk,i is the energy difference between the fc’th and 
the rth vibrational states of the same electronic state. 
The local energy separation 2tET^'^ for the n’th vibra¬ 
tional state is plotted in Eig. [^a). In Fig. [^a) we can 
also vaguely see frequency components corresponding to 


In Fig. [^b) we plot Ei - Wr for a number of energies Ei 
corresponding to vibrational states in different BO sur¬ 
faces. In practice the signal cos[(a;-Ai)T] does not exist 
for all T, and the peaks in Fig.j^b) are therefore not delta 
functions, but rather peaks smeared out around sets of 
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w’s and w^’s obeying the condition (42). If we focus on 
the absorption feature at ~ 11.5 eV we see from Fig. [^b) 
that the corresponding field-free state couples to at least 
three other vibrational states through the IR field. From 
this analysis we conclude that the interference mecha¬ 
nism resulting in T®/2 oscillations in the molecular ab¬ 
sorption lines of Fig. is very similar to the mechanism 
resulting in T^^^/2 oscillations in the fixed-nuclei ab¬ 
sorption lines discussed in Appendix]^ Other processes 
such as LIS and Autler-Townes splittings of the absorp¬ 
tion lines have not been found in the analysis of moving 
nuclei ATAS spectra of H 2 . 

We now investigate whether LIS at well-defined ener¬ 
gies are hidden in the rich molecular structure of ATAS 
spectra in H 2 , or a more fundamental nuclear motion 
effect suppresses such LIS in molecular systems. For 
this purpose we use a perturbative electronic three-level 
model similar to the one used in Eqs. (|T7|-(|^ for H 2 '*' 
with fixed nuclei, but now accounting for nuclear mo¬ 
tion. The model contains the electronic states (j)g{z;R) 
(ground state), R) (dark state) and (j)e{z-, R) (essen¬ 
tial, dipole allowed, excited state). The wave function of 
a general state in the model can be written as 

'i'{z,R,t) =Gg{R,t)4>g{z;R) + Gd{R,t)4>diz]R) 

+Ge{R,t)4'e{z]R) 

=Xg.o(^)<(>g(z;i?) 

+ ^ dAe,feCe,fe(t)Xe,fe(A)e“*^'''‘^*“^Ve(2:; R) , 

(43) 

where the nuclear wave packets Gg{R,t), Gd{R,t) and 
Ge{R,t) have been expanded in their vibrational eigen¬ 
states [see Eq. (30)]. In Eq. (43) the subscripts g,d 
and e label the three electronic states, while the sub¬ 
scripts k and I label the vibrational states corresponding 
to the electronic states e and d, respectively. In obtain¬ 


ing Eq. (43) we have approximated Gg{R,t) by its vi¬ 
brational ground state wave function Xg,o{R) use 
Ag o = 0. The approximation Gg{R,t) ^ Xg,o{R) is rea¬ 
sonable for the field parameters used in ATAS and is simi¬ 


lar to the approximation Cg{t) ^ 1 made in Sec. Ill A Un¬ 
der the approximation dfg{R) = df^{Ro) [see Eq. (34)] 


the expectation value of the dipole moment found using 
Eq. (43) is well approximated by [see Eq. (18)] 


{d-{t)) :^2Re 




dE, 


(Xg.oIXe.fc) Ce,fc(t)e 


e,k 


-iEc,k{t-r) 


(44) 


As in Sec. Ill A we approximate the NIR (or IR) field by 
a monochromatic field E{t) = E^ sm{ijjct) and evaluate 
the coefficients Ce^k{t) by first order perturbation theory 


„(i) 


it) = 


2i 


^dE, 


-f 

Ee,k ~ ^d,l 

- e~ 

~ I — U)p 


d,l 


(Xe,k\Xd,l)cd,l{A- 

(45) 


Inserting Eq. (45) into Eq. (44) shows that the ener- 


d,i 


gies Ee^k disappear in the final expression for the time- 
dependent dipole moment, leaving terms oscillating with 
frequencies Ed^i ±u!c- The distribution and population of 
vibrational states corresponding to the electronic state 
(j)e{z,R) does therefore not affect the dephasing of the 
dipole moment responsible for LIS. The distribution and 
population of vibrational states corresponding to the 
electronic state 4‘d{z,R) are, on the other hand, crucial 
for the existence of LIS as spectral structures at well- 
defined energies. The more dark vibrational states pop¬ 
ulated, the faster the resulting dipole moment will de¬ 
phase. Even in the situation where a finite number of 
discrete vibrational states are populated, and one would 
expect a partial revival of the time-dependent dipole mo¬ 
ment, the revival is often absent due to the finite dura¬ 
tion of the NIR pulse. We therefore expect LIS with 
well-defined energies to be suppressed in the molecular 
ATAS spectrum when more than a few vibrational dark 
states are populated. To test this hypothesis we return to 
H 2 ^ for a short notice. In H 2 ^ the BO curve correspond¬ 
ing to the dark electronic state investigated in Sec. Ill 
is dissociative, and a continuum of vibrational states are 
populated. We therefore expect the LIS corresponding to 
this dark electronic state to be suppressed in the ATAS 
spectrum. The ATAS spectrum S{uj,t) using the pulse 
parameters from Fig. [T] and a fixed delay of r = 0 is 
shown in Fig. [^a). We now construct an artificial H 2 ^ 
molecule by replacing the BO curve Ed{R) corresponding 
to the dark state with Eg{R- Rt) - [Eg{Ro) - Ed{Ro)], 
for Rt = 0 this BO curve is parallel to the ground state 
curve Eg{R). According to the Franck-Condon principle 
only the dark vibrational ground state will be populated 
when the two BO curves are parallel. We therefore expect 
LIS in the ATAS spectrum of the artificial H 2 ^ molecule 
with Rt = 0, shown in Fig. We see that this is 

indeed the case since spectral features of the expected 
width, an energy separation of two NIR photon energies, 
and centered around the energy of the dark vibrational 
ground state, appear in the spectrum. It should be noted 
that the coupling to the dark state curve has not been 
changed, and the results should therefore only be inter¬ 
preted qualitatively. In Figs. 0c) andj^d) the ATAS 
spectrum is shown for the artificial H 2 ^molecule with 
non-zero Rt- As Rt is increased, the LIS are damped. 
The reason for this damping is that the slope of the BO 
curve, corresponding to the dark state, evaluated at Rq, 
increases with Rt and more vibrational dark states are 
therefore populated for larger Rt- In H 2 the dark state 
curves of interest are far from parallel to the ground state 
curve, and hence we do not expect any significant LIS 
in the ATAS spectrum. In ATAS experiments using a 
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Figure 7. Normalized ATAS spectrum S{uj,r) [Eq. @] 
calculated at r = 0 for (a) H 2 ^ and the artificial H 2 ^ molecule 
(see text) with (b) Rt = 0, (c) Rt = 0.234 and (d) Rt = 0.625. 
The four insets show the BO curves Eg{R) and Ed{R)- In 
the insets Ro is indicated by a dashed gray line. The pulse 
parameters are as in Fig. but with a fixed delay r = 0 
between NIR and XUV pulses. 


dissociate on the excited dissociative BO curve. Tempo¬ 
ral information about this process can be obtained using 
only one pulse, i.e., without dressing of the system. In 
H 2 the broad absorption features of H 2 '^ were replaced by 
rich structures in the ATAS spectrum. For H 2 we focused 
on explaining the oscillations in the absorption lines, and 
found that these could be put into two different cate¬ 
gories; one category was oscillations on the timescale of 
half the IR period, the other category was oscillations on 
a much longer timescale. We argued that the oscillations 
on the timescale of half the IR period were very similar 
in nature to the oscillations seen in the absorption lines 
of the atomic-like H 2 ^ with fixed nuclei. We therefore 
conclude that these oscillations are due to interference 
between populated states of the same symmetry. The 
oscillations on longer timescales, however, are a purely 
molecular signature associated with transfer of popula¬ 
tion between vibrational states. Finally we found that 
the existence of LIS at well-defined energies in molecular 
ATAS spectra depends critically on the number of pop¬ 
ulated vibrational dark states. As a result clear, unam¬ 
biguous detection of LIS can only be expected in systems 
where the BO curves of the ground state and the involved 
dark state(s) are nearly parallel, or if one in another way 
can populate only a few vibrational dark states. 


molecular gas target, we therefore predict that LIS will 
be very difficult to detect. An exception to this general 
rule might be N 2 where the BO surfaces are remarkably 
parallel. It should be noted that Eq. (44) is valid only 
when the detuning is large, and first-order perturbation 
theory can be applied. The conclusions, however, remain 
the same if perturbation theory of a higher order is used. 
Further, our calculations show that there is no significant 
difference in the nature of LIS in systems with large and 
small detunings. 


VI. CONCLUSION 


The effect of nuclear motion in ATAS spectra has been 
elucidated. We performed calculations in ID H 2 ^ and 3D 
H 2 with fixed nuclei, and found that the spectra of these 
systems are very similar to atomic spectra. Moreover we 
investigated the same two systems including nuclear mo¬ 
tion. For the systems with moving nuclei, we saw that 
the nature of the BO surfaces, corresponding to excited 
electronic states, were crucial. In H 2 ^ all excited BO 
surfaces lead to dissociation, and as a result the ATAS 
spectrum shows very broad absorption features reflect¬ 
ing the damping of the time-dependent dipole moment 
by the dissociating nuclei. In the theoretical analysis of 
these broad absorption lines we found that, for specific 
field parameters, the width of the absorption lines were 
solely system dependent. The width of the lowest absorp¬ 
tion line in the ATAS spectrum could therefore be linked 
directly to the time it takes for a nuclear wave packet to 
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Appendix A: Interference between excited states in 
fixed nuclei H2^ 

In this appendix we model the interference process 
between excited states of the same symmetry (under 
0 -z) resulting in a T^^^/2 periodic modulation of the 

corresponding absorption lines in the ATAS spectrum. 
We focus on the interference fringes in the absorption 
line at 25.5 eV in Fig. corresponding to an energy 
difference Eq between the ground state and the fifth ex¬ 
cited state. At large negative delays the population of 
the fifth excited state (j)Q{z-,Ro) has two contributions 
for t > to- A direct contribution from the XUV pulse 
alone, and an indirect contribution from the XUV and 
NIR pulses acting in a combined manner. In the indirect 
process the XUV pulse populates the third excited state 
(j) 4 {z\ Rq) at time r, and at a later time to the NIR pulse 
excites population to the fifth excited state through the 
absorption of two photons. Keeping only terms that con¬ 
tribute to a dipole signal with frequency Eq we can write 
up the wave function of the system as 

'f{z;Ro,t) =0i(z;i?o) + 
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=(/)i(z;i?o) + (cd + Cide*"^) 

{t>to), (Al) 

where (j)i{z]Ro) is the ground state with energy Ei = 0, 
Ei,i = 4,6 are the energies of (I) 4 {z;Rq) and 4>^{z]R[)), Cd 
and Cid are the amplitudes corresponding to the direct 
and the indirect absorption processes, and the phase (p is 
given by 


p={EQ-E 4 )x{to-T). (A2) 

This phase is similar to the phase found in Ref. m where 
constructive interference conditions are investigated. To 
obtain Eq. we used that the population of the 

ground state is almost constant in time due to the low 
intensity of the XUV pulse. The time-dependent dipole 
moment in the simplified model is given by 


'O 

cos [Ee(t-T)] 

Cd cos [Ee{t - t)] + Cid cos [Ee{t - r) 




(t<T) 

(t < t < to) 
{t>to) 

(A3) 


In Fig. 1^ we compare the ATAS spectra constructed 
from dipole signals found from the full TDSE calculation 
[Eq. (113] and from Eq. (A3), respectively. Following 
Ref. [TT] we have chosen to to be located a quarter of 
a NIR period from the center of the NIR pulse. The 
amplitudes Cd and Cid can be found using perturbation 
theory, but the choice Cd = Cid = 0.5 suffices to illustrate 
the predictions of the model. As a result of the narrow 
spectral width of the NIR pulse, interference of the type 
seen in Eq. (A3) only occurs when the energy difference 


between two states of the same symmetry under z ->• -z is 
~ 2wnir. Consequently the associated oscillation periods 
are always ~ T^^^/2 with the period of a NIR pulse 
cycle. We note that the interference phenomena of the 
type discussed here, always affect both absorption lines 
involved (at E 4 and Eg), and a similar pattern is found 
in the absorption line corresponding to the energy E 4 . 
For this absorption line, however, the phase of Eq. ( |A2[ ) 
is {E4 - Eq) X (to - t) and the tilt of the fringes in the in¬ 
terference pattern will be opposite to that seen in Fig. [^ 
In Sec. jV] we see that interference phenomena similar to 
the ones investigated here are also present between vi¬ 
brational states in molecular systems, again leading to 
oscillations on the period of T^®/2. 
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Figure 8. (Color online) Normalized ATAS spectrum 
~S(u,t) [Eq. (§] calculated for the fixed nuclei model of 
H 2 ^. In the upper panel the Fourier transform of the time- 
dependent dipole moment is found from a full TDSE calcula¬ 
tion [Eq. ( |14[ )], while in the lower panel Eq. ( |A3| ) was used. 
The pulse parameters are as in Fig. [^ 
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